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1. Introduction 

As is well-known, the QCD pressure is an observable that plays a role in many contexts: 
besides being important for theoretical studies of the QCD phase transition and thermodynamics, it 
also has potential phenomenological relevance for cosmology and heavy ion collision experiments. 

At high temperatures, a useful approach for the computation of this observable is Dimensional 
Reduction [^, ^ ||]. It consists of replacing the full 4d theory with an effective 3d one including 
an adjoint Higgs field (Electrostatic QCD, "EQCD"). This theory can, in turn, be reduced to a 3d 
pure Yang-Mills theory (Magnetostatic QCD, "MQCD"). This strategy is useful first of all from the 
theoretical point of view, since it allows for a separation of contributions coming from the various 
scales that characterize QCD, namely T (hard modes), gT (soft modes) and g^T (ultrasoft modes). 
Moreover, it permits a study of the whole T-range of interest: the high-temperature region is usually 
investigated by means of perturbation theory while the low-temperature regime is explored via 
lattice simulations. There might, however, be a gap between the two regimes: on the perturbative 
side it is not possible to lower T too much because of the poor convergence [Q], while on the 
lattice side numerical limitations forbid simulations at temperatures higher than about STc |^]. 
Dimensional Reduction can overlap with both of these regimes and thus fill the possible gap. 

Within this framework, our first aim is to complete the determination of the order 0{g^) 
weak-coupling expansion of the QCD pressure: due to the presence of IR divergences [|6|], non- 
perturbative lattice measurements are needed at this order [Q], but their proper interpretation in 
the context of the full computation [§, ^ requires a conversion of the regularization scheme from 
lattice to MS. Second, the full Dimensional Reduction program requires the study of EQCD [|lO|], 
but the continuum extrapolations that enter at this stage turn out to be very delicate, and require the 
removal of lattice artifacts at 0{a) and 0{a^). 

Our aim is to compute these renormalization constants and improvement coefficients by means 
of Numerical Stochastic Perturbation Theory (NSPT), a procedure developed by the Parma group. 



2. NSPT basics 



NSPT has its origins in the concept of Stochastic Quantization [11] whose recipe is made up 
of two ingredients: the introduction of an extra coordinate, a stochastic time t, and an evolution 
equation of the Langevin type, 

^^-^.,(»), 

where T](x,f) is a Gaussian noise. Starting from this, the usual Feynman-Gibbs integration can be 
reproduced by averaging over the noise T], or more practically over the stochastic time t, that is 

Z-'l = limi |Jdf'(O[0,(x,f')]>, . (2.2) 

When dealing with SU (3) variables, the Langevin equation needs to be modified into 

^,U^ = -i(wS[Ur^] + ^)u^, (2.3) 

in order to assure the correct evolution of the variables within the group. 
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In this framework, perturbation theory can be introduced by means of the expansion [12] 



(2.4) 



where go is the bare gauge coupling. This gives a system of coupled differential equations that 
can be solved numerically via a discretization of the stochastic time t = nT, where T is a time step. 
In practice, we let the system evolve according to the Langevin equation for different values of T, 
average over each thermalized signal (this is the meaning of the above-mentioned limit f — > oo), and 
then extrapolate in order to get the T = value of the desired observable. This procedure is then 
repeated for different values of the various parameters appearing in the action. 



3. Renormalization of the Magnetostatic sector: setup 

The MQCD contribution to the pressure can be written as [| 



6 "A 
SM 



{Ak) 



43 157 2\, 
12 - 768^ '1"- 



+ Bg + 0(s] 



(3.1) 



with Nc the number of colours, cIa = — I, At the MS scheme scale parameter, gM the gauge 
coupling and Be the constant we ultimately want to determine. It can be shown that for Nc = 3 [pj|]. 



5g= 10.7 ±0.4 -5l(1)+Sms(1) + 



43 157 



1 



+ ln2 + 21nA/c 



(3.2) 



12 768 J \3 

where the numerical value is non-perturbative and follows from lattice simulations |^] (for the 
corresponding numbers atA^c ^ 3, see ref. Q), while Bms(1) = -2.16562591949800919016 has 
been determined as a result of extensive continuum computations [15, |T^]. The remaining un- 
known, Bl(1)> can be expressed as [ pj| ] 



lim 6o' 



1 Tr[Pi2l 



up to 4-loop 



_£2_ C3_ C4 J_ 



(3.3) 



where the argument "1" corresponds to Feynman gauge, jSq = 2Nc/aglj,a is the lattice spacing, 
Tr[Pi2] is the trace of the elementary plaquette in the 1-2 plane, and m is a gluon mass the has been 
introduced as an intermediate IR regulator. The coefficients ci , . . . , C4 are all known [ 17 , 1^, |19[ 0] . 



To evaluate eq. (p.3|), gauge fixing and mass terms need to be introduced: 

Z = y [D0] exp {-S„-SaF- S^p) , 
where we assume the use of lattice units (i.e. a = 1), and 



4N, 



^ x,A " II 



ln(-£5^D^[0]+m^ 



(3.4) 

(3.5) 

(3.6) 
(3.7) 
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Here we have followed the conventions of ref. pO|], writing in particular = exp(/^^), (^^ = 
<pflT^, with the normalization Tr[T^T^] = 5^^ /2. Moreover m is the common gluon and ghost 



mass and is the discrete Faddeev-Popov operator, given by [ |20| ] 



1 _|_ i<T>„ L<T>2 

12^ 



1 

720^ 



1 



30240 



06^ + 0(0^) 



(3.8) 



with = (^t^F^, where [F'^Jgc = —if^^^ are the generators of the adjoint representation. Details 



about the treatment of the Faddeev-Popov determinant can be found in ref. [ ]13p . 

The procedure then consists of measuring the plaquette for different lattice sizes at fixed mass, 
then extrapolating towards infinite volume, and repeating for different values of the mass. Finally, 



after subtracting the logarithmic divergence, the zero-mass extrapolation of eq. (3.3) will provide 
the quantity we want to measure. It is important to perform first the extrapolation in volume and 
then in mass, because the opposite order would result in having the finite size as the IR regulator 
and not the mass as desired. 



4. Renormalization of the Magnetostatic sector: results 

As just stated, the first step is the extrapolation in volume. While the analytic behavior is 
known at 1-loop level (and this has given us a useful crosscheck), this is not true at 4-loop level, so 
that we have to rely on effective fits. Some of them are shown in Fig. 1. 

In order to be as conservative as possible, we opted for fitting a constant to those points that 
do not seem to show any volume dependence within errorbars. To check whether this approach is 
reliable, we employed it for the first three loops, and then performed the zero-mass extrapolation to 
see whether the already known coefficients c\, C2 and C3 are recovered. Table 1 confirms that this 
indeed is the case. 



Coefficient 


Our extrapolation 


Known result 


Cl 


2.672(8) 


2.667 


C2 


1.955(16) 


1.951 


C3 


6.83(10) 


6.86 



Table 1: Comparison of our extrapolations with the known results [|l7], |9[ ^ for the first three loops. 



The same procedure was subsequently applied at the 4-loop order where, however, the IR 
divergence needs to be subtracted before taking the zero-mass limit (see Fig. 2). We then performed 
polynomial extrapolations involving a different number of points and degrees of freedom: the final 
results we get for 6^(1) is [ p^ 

fiL(l) = 13.8 ±0.4 , (4.1) 



which gives, once inserted into eq. ( |3.2D , 

Bg = -0.2 ± 0.4^^^) ± 0.4(^^^^) . (4.2) 
Here "MC" labels the result of Monte Carlo simulations |^. 
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Figure 1: 4-loop plaquette vs mL at fixed mass 
(m — 0.2): the solid and the dashed lines cor- 
respond to different combinations of a constant 
plus a sum of negative exponentials and negative 
powers of mL; the horizontal line fits a constant. 



Figure 2: Infinite-volume 4-loop plaquette vs m: 
circles still contain ln(m) while diamonds are ob- 



tained after its subtraction (cf. eq. (3.3)) 



5. Renormalization and improvement in the Electrostatic sector: setup 

The EQCD action in the continuum is given by 

Se = j d^'^|iTr[F,2(x)] +Tr[A-,AoW]' + m2Tr[A2(^)] + A(Tr[A2(^)])'| , (5.1) 

where Fij{x) is the 3d field strength tensor, D,- the covariant derivative and Aq = Afr^, with 
normalised as before. Once again the strategy is to obtain the MS scheme result starting from 
lattice measurements. Apart from the plaquette expectation value, this theory has however more 
condensates that play a role. In particular, derivatives with respect to the dimensionless variables y 
and X, defined as j = /g\, x = X /g\ (with gE the EQCD gauge coupling), produce condensates 
quadratic and quartic in Aq [1Io|]. Subtracting the proper counterterms [||], and 0{ a) or 0{a^) 
effects, which become important in the range of large y where connection to the weak-couphng 
expansion can be made, we can extrapolate to the continuum, and finally obtain the pressure by 



integration []10|]. 

The EQCD lattice action is given by 



Sin 



PL 



•X, i<j 



1 



Re Tr [Pij {x)])-2 1,, Tr [0 (x) Ui (x) (x + /) uj (x)] + 



+ iJ a(j8,A,jiatt)Tr[(/.2(x)] +x(tx[^\x 



(5.2) 



where now j8 = 2Nc/ag\, Ui is the link variable, = Aoy67j8, and |21] 
a(/3, A, jMt) = 6| 1 + Ijiatt - (6 + \X^)'-n^^^^- 



3 



(10Aj8 - |A2/32) (In j8 + 0.08849) + M^AjS + 36.130 



(5.3) 



5 



Renormalization of infrared contributionsto the QCD pressure 



C. Torrero 



The first condensate we want to measure is tlie derivative witli respect of yx^ix'- apart from a 
rescaling factor, it is equal to (Tr[AQ]) whose lattice counterpart has a perturbative expansion given 
by 

These include the counterterms and lattice artifacts mentioned above: some of the coefficients (cfoo> 
diQ, d\\, d2i, d22) have already been estimated while the other ones (especially J20 and J30) are 
what we aim at computing by means of NSPT. We will again measure the observable for different 
values of the lattice extent L and the parameters InjS and yMt, and carry out an extrapolation in L 
(at fixed lnj8 and yun) to get the infinite- volume results, from which we infer the behavior of the 
coefficients when varying the other variables. 



6. Renormalization and improvement in the Electrostatic sector: first tests 

The statistics we have collected so far is sufficient just to check the reliability of this approach. 
A first test is to compare our numerical estimates for the coefficient doo at fixed L and y^t with the 
analytical results: this comparison is shown both in Table 2 and in Fig. 3 and appears satisfactory. 



L 


Exact result 


NSPT estimate 


5 


0.6861 


0.6867(11) 


6 


0.6833 


0.6845(8) 


7 


0.6825 


0.6837(6) 


8 


0.6822 


0.6825(5) 


9 


0.6821 


0.6822(4) 


10 


0.6821 


0.6829(4) 



Table 2: Comparison between exact and NSPT values of doo, for y\att = 1 0. 



As a second check, we inspect how our finite-volume data approach the infinite-volume limit 
at those orders for which we have a direct "exact" estimate of this limit: an example is given in 
Fig. 4 for the coefficient d2i (whose infinite-volume value is 1.4072). Once again, the behavior 
looks encouraging; the same is observed for the terms not shown here. 



7. Conclusions and prospects 

While our determination of the renormalization constant related to the 0{g^) contribution to 
the QCD pressure from the Magnetostatic sector has recently been completed [|T3|], there is still 
work to do as regards the Electrostatic contributions. It is important to finalise this task, since the 
EQCD result has a wider range of applicability than the MQCD result alone. The first tests have 
produced encouraging results, so that there is every reason to believe that the determination of the 
most important new coefficients (d20 and ^30) is also feasible, at least in a certain range of jiatt. 



With these results, the program initiated in ref. [10] could finally be carried out to completion 
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Figure 3: The coefficient c/qo vs L; exact (dia- 
monds) and NSPT (circles) results. 
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Figure 4: The coefficient d2\ vs L. The expected 
infinite-volume value is 1.4072. 
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